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We discuss recent work on the static and dynamical properties of the asymmetric exclusion process, 
generalized to include the effect of disorder. We study in turn: random disorder in the properties of 
particles; disorder in the spatial distribution of transition rates, both with a single easy direction and 
with random reversals of the easy direction; dynamical disorder, where particles move in a disordered 
landscape which itself evolves in time. In every case, the system exhibits phase separation; in some 
cases, it is of an unusual sort. The time-dependent properties of density fluctuations are in accord 
with the kinematic wave criterion that the dynamical universality class is unaffected by disorder if 
the kinematic wave velocity is nonzero. 
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I. INTRODUCTION 

Simple models of particles moving on a lattice have 
shed a good deal of light on the general question of collec- 
tive effects in transport. Even the simplest sorts of inter- 
actions between particles can induce collective phenom- 
ena involving large numbers of particles, ranging from 
the motion of density fluctuations as a kinematic wave, to 
the occurrence of phase transitions to states with macro- 
scopically inhomogeneous density (jams). Models which 
display these phenomena include the asymmetric simple 
exclusion process (ASEP), the simplest model which in- 
corporates directed motion and mutual exclusion Q,Q,0|, 
and its generalizations which include more realistic fea- 
tures which mimic vehicular traffic 0. 

Here we are concerned with the effects of extensive dis- 
order on interacting-particle transport in one dimension. 
As in equilibrium systems, disorder has a strong effect on 
the properties of driven systems which reach a nonequi- 
librium steady state. However, in contrast to equilibrium 
systems, relatively little is known, in a general way, about 
the effects of quenched disorder on nonequilibrium sys- 
tems, as the absence of detailed balance, together with 
the breaking of translational invariance, makes even the 
determination of the steady state weights difficult in gen- 
eral. In this backdrop, the study of simple models of dis- 
ordered nonequilibrium systems has proved valuable, as 
it reveals how interactions, drive and disorder combine to 
produce new types of states; the simplicity of the models 
allows for a detailed characterization and understanding. 

In this paper, we will review some results obtained for 
a disordered system of particles moving stochastically on 
a one-dimensional lattice. We restrict ourselves to the 
disordered ASEP, in which interactions between parti- 
cles enter through rules for the hopping rates - the hard- 
core constraint is modelled by forbidding any move which 
would lead to more than one particle per site. We will 
consider three types of disorder. The first is particle-wise 
quenched disorder, in which different particles have dif- 
ferent properties which do not change in time. Secondly, 
we will consider space-wise quenched disorder, in which 
particles do not differ from each other, but are subject to 
hopping rates which are randomly distributed in space. 



Finally, we will turn to dynamical disorder where parti- 
cles move in a disorder-induced landscape, which is itself 
evolving in time. Interestingly, there are inter-relations 
between the effects of the three types of disorder, which 
we discuss below. 

It is clear that extensive disorder at the microscopic 
scale would induce inhomogeneities in the particle den- 
sity on a similar scale. Furthermore, disorder can also 
induce variations of the density on a macroscopic scale, 
associated with phase separation p| . The models consid- 
ered here share this feature, but in some cases the phase 
separated state has quite unusual properties. We also 
discuss the time-dependent properties of these systems, 
both in the steady state and while approaching it. 



II. PARTICLE- WISE DISORDER 

Consider a system of cars on a single-lane road, too 
narrow for overtaking to be possible. Each car has an 
intrinsic maximum speed, which however may not be 
achieved due to a slowly moving car in front of it. It 
is then evident that the slowest car in the system will be 
trailed by a number of intrinsically faster cars, and will 
see a larger than average headway in front of it. The 
question arises: In the limit of an infinitely long road 
with a finite density of cars, will this headway be finite 
or infinite? Is there a phase transition between these two 
possible types of behaviour as the density is changed? 

This question motivates the study of a collection of 
particles with a different intrinsic rate of motion for each. 
In a simple lattice model, studied in and Q, N par- 
ticles reside on the sites of a one-dimensional ring with 
L sites. The dynamics involves random sequential up- 
dating of configurations. An elementary move consists of 
an attempted rightward hop of a particle, say the fcth, 
with a particle-dependent hopping rate u(k); the move is 
actually implemented only if the site immediately to the 
right is unoccupied. A configuration is specified by the 
set of particle locations {yk}- One may instead consider 
the set {mfc}, where rrik — J/fc+i — 2/fc — 1 is the head- 
way or gap between the fcth particle and the one ahead 
of it. On a periodic ring, a typical set {rrik} represents 
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N equivalent configurations which are translational shifts 
of each other. Evidently the overall density p = N/L is 
conserved by the dynamics. 

Disorder enters through the selection of hopping rates 
u(k), each of which is drawn independently from a dis- 
tribution Prob(u). Especially interesting is the case with 
a power law variation near a minimum cutoff speed: 
Prob(u) = [(n + 1)/(1 - u ) n+1 ](u - u Q ) n with u in the 
range [u , 1]. Here n is an index which characterizes the 
decay of the distribution function near the cutoff. As we 
shall see below, this system shows a phase transition as 
the density is decreased, from a state with finite head- 
ways ahead of each particle, to one in which the headway 
ahead of the slowest particle is infinitely long, or more ac- 
curately, a finite fraction of the of the lattice size L. 



A. Mapping to Zero Range Process 

The analysis of this system is aided by exactly map- 
ping this system with particle-wise disorder to another 
system with space- wise disorder 0, Q ■ View the particle 
index k as labelling urns arranged in a one-dimensional 
sequence (Fig. 1). The occupancy of urn k is taken to 
be the length of headway mj. in front of the fcth particle 
in the original particle problem. Now, a single rightward 
hop of a particle reduces the headway in front of it, and 
augments it behind. In the urn representation this cor- 
responds to an elementary move of a particle from urn k 
to urn k — 1 (Fig. 1). 



Here Z is the grand partition function, rn^ is the mass 
at site k, and z is the fugacity which is determined 
by requiring J2k TOfe ~ M. For a given realization 
of disorder {u(k)}, the density p = M/N is given by 
TV -1 J2k z l (u(k) — z). As pis increased, the model under- 
goes a phase transition We use the self- averaging 
property to replace the sum over k by an average over 
Prob(u), separating out the lowest u: 



P : 



du 



/(«) 



(2) 



where u m i n is the lowest realized value of u{k)\ u m i n 
approaches u in the thermodynamic limit. The fugacity 
z increases with p until one reaches a critical density 
Pc — Uo(n + l)/ n (l — u o) beyond which z gets pinned to 
u a . Beyond this point, the system enters a high-density 
phase in which the excess mass M — p c N is 'condensed' 
at the site with the lowest hopping rate in the ZRP. The 
transition is analogous to Bose-Einstein condensation in 
the ideal Bose gas or a system of bosons in a random 
potential [Icj . 

In the particle-wise disordered model, the condensate 
represents an infinite headway in front of the slowest car, 
implying an infinitely long jam behind it. Such a jam 
implies spatial phase separation, and exists only for low 
enough density; at higher densities, blocking effects be- 
tween cars reduce the contrast in speeds to such an extent 
that headways are interspersed throughout the system 
and there is no phase separation. 
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FIG. 1: Particle- wise disordered model and the Zero-range 
process 

This model is a special case of the zero range process 
(ZRP) which describes the motion of particles between 
wells or urns, with (in general), occupancy dependent, 
site-dependent hopping rates u(k,mf.) PUg. The ZRP 
has the virtue that the steady state measure P S {C) is 
known. In our case, where the rates depend on sites but 
not on occupancies, the steady state weights are partic- 
ularly simple. In a grand canonical formulation of the 
problem we have 



P S (C) 



z 



N 

n 

fc=i 



u(k) 



(1) 



B. Time-dependent Properties 

It is interesting to ask about the dynamical properties 
of the system, both in the steady state, and also describ- 
ing relaxation to the steady state. In steady state, it 
is evident that the mean current between every pair of 
neighbouring sites in the ZRP is the same. It is less ob- 
vious that the full distribution of time intervals between 
particle transfer events is also identical on each site. This 
property follows from the application of Burke's theorem 
from queueing theory [TH Il2j | to the problem at hand. 
This theorem states that in a reversible birth-death pro- 
cess, a Poisson input of births results in an identical Pois- 
son distribution of deaths. Identifying the occupancy of 
each urn of the ZRP with a population undergoing such 
a process then implies the result, as the output from one 
urn constitutes the input for the next. In the original par- 
ticle problem, this means that the time history of every 
particle follows an identical Poisson distribution despite 
the fact that each particle has different intrinsic speeds. 
This holds in the thermodynamic limit in both the dis- 
ordered and the phase separated state 01 ■ 

Next, let us ask about the behaviour of statistical fluc- 
tuations of the density in the steady state. In a homo- 
geneous driven system with a density-dependent current 
J(p), density fluctuations are transported around the sys- 
tem as a kinematic wave with speed c = dJ /dp, in time 
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N/d 0,113. This wave of fluctuations decays on a time 
scale ~ N z (z > 1), implying that it circulates several 
times around a system with periodic boundary condi- 
tions before it dissipates. In the disordered ZRP, there 
are different sorts of behaviour in different phases. For 
p < p c and for p = p Cl n > 1, the wave speed c is 
nonzero. There is clear evidence of the wave in simula- 
tion results for the fluctuations in the integrated num- 
ber of particles which move past a site. This correlation 
function shows pronounced oscillations with period N/c, 
which ultimately damp down |l5j| . 

In [3 it has been argued that quenched disorder does 
not have a significant effect on the damping of fluctu- 
ations in driven systems if the kinematic wave speed c 
is nonzero. This 'kinematic wave criterion' implies that 
there is then no change in the universality class of the 
decay of density fluctuations, which should continue to 
follow the behaviour in the absence of disorder. The rea- 
son is that in such a case, each density fluctuation does 
not stay long in the vicinity of a local patch of disor- 
der. To that extent, the effect of quenched disorder on a 
ballistically moving density fluctuation of spatial extent 
Air is akin to temporal noise, and we should expect the 
lifetime of the density fluctuation would continue to be 
~ (Ax) 2 , where z is the dynamical exponent in the ab- 
sence of disorder. Monte Carlo simulations [15( confirm 
that when c ^ 0, z equals 3/2, the value for a driven sys- 
tem in the absence of disorder [l^. At p — p c , the speed 
c vanishes for n < 1. In this case, disorder is expected to 
have a strong effect, and perhaps change the universality 
class away from that in its absence. The determination 
of this class remains an open problem. 

For p > p c , the system has density p c in the bulk, 
with the surplus mass (M — pN) being on the slowest 
site. Thus, density fluctuations would be expected to 
move as a kinematic wave in the bulk so long as n > 1. 
However, this wave would be absorbed at the condensate 
site and not be able to circulate, implying that correlation 
functions in this phase would be non-oscillatory. 

Now let us turn to the kinetics of approach to the 
steady state, starting from a state with uniform density. 
In the disordered phase (p < p c ), the kinetics is gov- 
erned by the same dynamical exponent z that operates 
in steady state, and the relaxation time is of order L z , 
with z = 3/2. In the condensate phase, however, a dif- 
ferent physical process, namely coarsening, governs the 
approach. Initially, particles hop out of relatively fast 
sites quite quickly, and get trapped at slow sites in the 
immediate neighbourhood. At moderate times, there is 
thus a finite density of noticeably large aggregates at such 
slow sites, which then relax by releasing their mass to yet 
slower sites ^tJ- Thus the mass at all such slow sites, ex- 
cept for the very slowest, shows a nonmonotonic variation 
in time. The mass at the slowest (condensate) site grows 
as Mi ~ tP, with = (n + l)/(n + 2), a result first de- 
rived for a model with deterministic dynamics [l^, but 
which remains true even in the stochastic model under 
consideration Hi Kl [la . This translates into a statement 



about the growth of headway lengths £ in the particle 
model 

e(t)~t&*. (3) 

A similar coarsening description for headways is also 
found to hold in a lattice gas model of traffic, which in- 
corporates features of acceleration and slowing down, and 
in which different drivers have different propensities for 
random braking |20j . 

Furthermore, the problem has also been studied us- 
ing open boundary conditions, with injection and ejec- 
tion rates specified |2l|. In the absence of disorder, it is 
known that the open system exhibits phase transitions 
between a low-density phase, a high-density phase, and 
a high-current phase, as the injection and ejection rates 
are varied [2^,|2j|. With particle- wise disorder, there is 
a shift of the phase boundaries owing to the breaking of 
particle-hole symmetry. Within the low-density phase, 
there is a regime in which the drift velocity of all par- 
ticles is determined by the slowest [2l|, as with perodic 
boundary conditions. 

Finally, we discuss some recent work on two-way traf- 
fic |2(| |23 ■ Consider particle- wise disorder in which the 
majority of particles move preferentially rightward, while 
a minority (a fraction /) move peferentially leftward. For 
each particle, the ratio of the larger hopping rate to the 
smaller one, which is a measure of bias, is taken to be 
fixed, but the direction of hopping is random. The prob- 
lem was mapped onto a ZRP, which inherits the random 
easy-direction property. This system has much in com- 
mon with the spatially disordered ASEP with bidirec- 
tional bonds, to be discussed in the next section. In par- 
ticular, as in that case, the current decays as an inverse 
power of the size, with the power depending continuously 
on both / and the bias. 

III. SPACE-WISE DISORDER 

Consider the transport of interacting particles which 
are driven through a randomly disordered one- 
dimensional path. To address the question of what sort of 
macroscopic states result, we restrict ourselves to the dis- 
ordered asymmetric simple exclusion process. Disorder 
enters through the transition rates, which are assigned 
randomly to the bonds between neighbouring sites, e.g. 

describes the rate of attempted hopping from site 
i to i + 1, while ttj+i,i is the attempt rate in the reverse 
direction. A hop to a site actually occurs only if the site 
in question is unoccupied. Results obtained for a par- 
ticular disorder realization R = {ui.i+i, iti+i.i} are then 
to be averaged over R using a pre-specified probability 
distribution, taken to be the product of identical and in- 
dependent distributions across all bonds. The larger of 
the two rates on each bond defines the easy direction on 
that bond. Below we will distinguish between two situa- 
tions: unidirectional, in which the easy direction of every 
bond is the same, but the strength is a random variable 
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((a) and (b) in Fig. 2) and bidirectional, in which the 
easy direction is itself a random variable ((c) and (d) in 
Fig. 2). 

In the unidirectional case, we will further specialize to 
only the case of no backward hopping (u^+i.i = 0): 

Prob(u hl+1 ) = (l-f)S(u i s +1 -u 1 ) + f5(u i s +1 -u 2 ) (4) 

with U2 < u\. The fraction / and relative strength r = 
U2/111 of weak bonds specify the extent of disorder. 

In the bidirectional case, the easy direction of each 
bond is taken to vary randomly but the relative hopping 
rate A against and along the easy direction on each bond 
is taken to be the same. 



(a) Unidirectional disorder (Low density) 



Prob(b ht+1 = ^±M) 

Ui,i+1 



(l-f)6(b i , i+1 -X)+f5(b i , i+1 -\- 1 

(5) 

Here / is the fraction of backward pointing bonds. 

The large scale behaviours that result depend on the 
type of disorder as well as on the overall density p. Figure 
2 shows a schematic depiction of a segment of the disor- 
dered lattice, the corresponding variation of the current 
J with system size L, and the schematic spatial variation 
of the mean occupancies pi = (rii) in a particular real- 
ization of disorder in different cases. These are discussed 
below. 

(a) With unidirectional disorder and a particle density 
that is sufficiently far from 1/2, i.e. \p — 1/2| > A, the 
density is homogeneous on the macroscopic scale. The 
value of A depends on the concentration / and relative 
strength r = u 2 /ui of the weak bonds. The current J 
approaches a nonzero value in the thermodynamic limit. 

(b) With unidirectional disorder and p close to 1/2, 
i.e. \p — I/2| < A, the state exhibits phase separation; 
the system is characterized by two distinct values of the 
density, each extending over a macroscopic region. As in 
(a), J has a finite value in the thermodynamic limit. This 
value of the current is the same for all p in the regime 
|p-l/2| <A. _ 

(c) With bidirectional disorder, but with an asym- 
metry in the number of forward and backward bonds 
(/ 7^ 1/2), the potential has local minima and an overall 
tilt. On a macroscopic scale, the system has two distinct 
values of the density close to 1 in one region and close to 
in the other. The current J falls with increasing system 
size L as an inverse power of L. 

(d) With bidirectional disorder and equal numbers of 
forward and backward bonds (/ = 1/2), the potential has 
local minima, but no overall tilt. Consequently the mean 
current is zero. The density profile shows large sample- 
to-sample variations, but every sample has a high-density 
region that extends over a finite fraction of the system 
size. 

Let us turn to an explanation of the features described 
above. Central to the discussion is the observation that 
while there are strong spatial variations in the local den- 
sity (rii) = p i} the steady state value of the current J is 
uniform across all bonds. 



(b) Unidirectional disorder (Intermediate density) 




(c) Bidirectional disorder (Tilted) 



V 




(d) Bidirectional disorder (Untilted) 




L i 

FIG. 2: Space-wise disorder: depiction of possible behaviours 

A. Unidirectional disorder 

The time averaged current in bond (i, i + 1) is J^j+i = 
Ui t i+i(rii(l — m+i)). Since the exact steady state weights 
of particle configurations {rii} are not known, one cannot 
compute the two-point correlation function involved. To 
proceed, we employ a mean field approximation which 
ignores correlations betwen site densities in the steady 
state, and replaces {nifij) by (rii)(nj) = ptpj, but does 
allow for spatial variation of pi, which is crucial. To 
find the mean field solution {pi} we may proceed in two 
ways, (i) For a given value of J, iterate the set of equa- 
tions pi + i = 1 — J/v,ij+\pi around the periodic chain, 



5 




P 



FIG. 3: Current-density plots 



until convergence is achieved, (ii) For a given value of 
the overall density, evolve local densities (in fictitious 
time) through pi(t + 1) = p t (t) + ji-i,i(t) - ji,i+i(t) 
where the instantaneous currents are given by ji t i+i(t) = 
u i,i+iPi{t)[l — Pi+i(t)]- The densities then converge to 
time-invariant values, and the corresponding current is 
uniform across bonds. As discussed in |24|, method (ii) is 
preferable, as values of J above a threshold value are not 
allowed, and method (i) does not converge well near the 
threshold. The resulting mean field density profile can 
be tested against the results of direct Monte Carlo sim- 
ulation. The agreement is remarkably good, in that the 
mean field approximation tracks every one of a macro- 
scopic number of 'minishocks' and correctly predicts their 
positions, though it does not accurately reproduce their 
shapes [24j . 

The principal result is the existence of a density 
range A in which the current has a constant, density- 
independent value (regime (b)), and where the system 
exhibits phase separation of densities [24| . A simple qual- 
itative understanding of both these features can be ob- 
tained by referring to Fig. 3. The argument relies on 
a conjectured maximum current principle, which states 
that for a given value of the overall density, the system 
settles into a state which maximizes the mean current. 
In our system, let us suppose that the density in each 
stretch of like bonds is uniform, and label by 1 and 2 
the stretches of strong and weak bonds respectively. The 
two parabolas in Fig. 3 are the corresponding J versus 
p curves for the pure 1 and pure 2 systems. Since the 
current is uniform in the disordered system, a value of 
current such as shown by by the upper dashed line e/ in 
the figure is ruled out, as such a value cannot be sustained 
in a 2-stretch. By contrast, a value of J corresponding to 
say ad is quite possible, as the (1,2) densities can assume 
the values (p a ,Pb) if the density is relatively low, or al- 
ternatively (p c , pd) if the density is relatively high. Such 
choices correspond to the regime \p— 1/2 1 > A, and result 
in a density profile that is homogeneous on a macroscopic 
scale as strong and weak bonds are interspersed at the 
microscopic level. If the density is too close to half-filling 
\p— 1/2 1 < A, the system adjusts by keeping the current 



constant at its maximum allowed value, and replacing 
{PaiPb) m (1)2) stretches by (p C iPd)i m a finite fraction 
of the full lattice. The resulting state then shows phase 
separation of the density. The location of the high den- 
sity region is decided by the longest stretch of weak bonds 
in the system, as this stretch (whose length is of order 
InL) mimics the effect of a single defect weak bond in an 
otherwise pure system, which is known to induce phase 
separation in a finite range of density [2^, [2j| . 

Interestingly, it is possible to obtain upper and lower 
bounds |f| on the value of the threshold density for phase 
separation p c = 1/2 — A in terms of the ratio of bond 
strengths r = u%/ux < 1 and the fraction / of weak 
bonds. An upper bound is obtained on noting that the 
maximum current that can flow through a stretch of slow 
bonds decreases as the stretch length increases. It is then 
plausible that the current would be the least in the fully 
segregated limit where all slow bonds form a single large 
stretch |24| . In such a system consisting of two connected 
homogeneous parts, it is easy to figure out the way a 
given overall density would distribute itself over the two 
portions. A lower bound on p c was obtained [j| by argu- 
ing that the current in the disordered ASEP is bounded 
above by the current in a corresponding ZRP with the 
same disorder distribution and an equal number of par- 
ticles. 

The results discussed above hold for the case of binary 
disorder in the transition rates. Similar phenomena, e.g. 
the plateau in the J — p curve over a finite range A of 
density, have been found also for more general, contin- 
uous distributions of the transition rates pffll l3ll | . The 
problem has also been studied using open boundary con- 
ditions, specifying injection and ejection rates. An ar- 
gument based on moving shock fronts [s3 suggests that 
in the presence of disorder, the high-current phase ex- 
tends over a larger region of the phase diagram [3lJ than 
in the pure system [22j. Further, a numerical study of 
the first order transition between the low and high den- 
sity phases shows that the transition remains first order, 
but the location of the transition shows strong sample to 
sample fluctuations which do not seem to damp down in 
the thermodynamic limit |33| . 

Going beyond the ASEP, the effect of slow sites on a 
lattice model of traffic was studied by numerical simu- 
lation in j^]. Interestingly, a regime with two current 
plateau regions was observed, and the system was ob- 
served to switch between the two branches. 

To conclude this section, we mention an intriguing 
symmetry that holds for the disordered ASEP when no 
backward hopping is allowed. A typical configuration 
of quenched disorder is not symmetric under space in- 
version, and consequently nor is the steady state density 
profile. In spite of this, the magnitude of the steady state 
current was observed to be invariant under reflection |24j . 
Subsequently, this invariance property has been proved 
for both periodic [3{|,|3(j and open [3|| systems. 
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B. Bidirectional disorder: With tilt 

When the easy direction on a bond is a random vari- 
able, one needs to distinguish between the cases / < 1/2 
and / = 1/2, corresponding to a tilted and untilted po- 
tential respectively (figs. 2(c) and (d)). We may imag- 
ine assigning arrows to denote the easy direction on the 
bonds of a periodic ring; a fraction (1 — /) point right- 
ward, while a fraction / point leftward. A particle-hole 
exchange occurs in the direction of the arrow with rate 
u, and in the opposite direction with rate it A with A < 1. 

In the tilted case, there is an overall tendency for the 
particles to move rightward, but the question is whether 
there is a nonzero current in the thermodynamic limit. 
Local barriers to the rightward motion occur in the form 
of backbends which consist of a chance agglomeration of 
successive left-pointing arrows. An upper bound J\ on 
the current that can be carried by a backbend of length 
I is obtained by considering the boundary conditions, 
p = 1 at the leftmost end of the backbend and p = 
at the rightmost end, which force the largest possible 
current through it Within a mean field approximation, 
the density profile, and hence the current, can be cal- 
culated |33 with the result J/ ~ A The factor ~ is 
a consequence of particle-hole symmetry which implies 
that the particle density extends up to half-way through 
the stretch, so that the topmost particle needs to be acti- 
vated across hi sites. An exact calculation [2!| confirms 
the exponential decay of current with I. The current in 
the full system of size L is limited by the length l*(L) 
of the longest backbend. Since the probability of the oc- 
currence of I successive backward bonds is / , we may 
estimate I* from Lf l ~ 0(1). The estimate of the cur- 
rent A' I" 1 can be written as 

J{L) ~ L- 9 ' 2 (6) 

where the exponent 

9 = InA/ln/ (7) 

depends explicitly on the bias and fraction of reverse 
bonds. This dependence has been confirmed by Monte 
Carlo simulations [24| . 



Since each local slope of the height is equal in magnitude 
but random in sign, the full height profile is isomorphic 
to the trail of a random walker, with the height being 
the displacement of the walker and i being time. Thus 
the problem reduces to that of hard-core particles in a 
Sinai potential. In the limit of large bias (T — + 0), the 
system approaches the ground state, which can be found 
as follows. Make a constant-height cut of the height pro- 
file, fill particles up to that point, and then raise the 
level of the cut step by step and fill in particles till 
they are exhausted. The lengths £ of particle clusters 
(which are present below the topmost filled level) arc 
then determined by returns to the origin of the walker, 
and their distribution follows an asymptotic power law 
Prob(£) ~ £~ z / 2 . This distribution has an infinite mean, 
and in a system of size L an estimate based on extremal 
statistics shows that the largest cluster size is gL. The 
prefactor g shows strong fluctuations from one sample 
to another, but is always of order unity. Thus a typical 
configuration consists of at least one 0(L)-sized cluster 
in addition to ~ \[L clusters of size ~ \J~L. 

In j33, a similar model (the Coarse-grained Depth or 
CD model) was analysed, and the disorder-averaged two- 
point correlation function Cl(x) = {nini-\. x )i, — p 2 was 
computed, using the cluster size distribution Prob(£) as 
input. The result is that Cl{x) is a scaling function of 
argument x/L: 

C L {x) = Y{x/L) w m 2 [l-a(x/L) a ], x/L«l (10) 

with intercept m a = 1 and a = 1/2. 

The argument of the scaling function indicates that 
it describes a phase-separated state 39]. The finite in- 
tercept m 2 is a measure of long-range order. The cusp 
singularity is a manifestation of ill-defined interfaces, as 
it is known that a phase-separated system with well de- 
fined interfaces has a = 1 , in concordance with the Porod 
law [3^|. We conclude that the state exhibits an unusual 
sort of phase separation, with unusual interfacial charac- 
teristics and strong sample to sample fluctuations. 

D. Time-dependent properties 



C. Bidirectional disorder: Untilted 

When there is no overall tilt, the mean current is iden- 
tically zero. The problem then reduces to the equilib- 
rium problem of hard core particles at temperature T 
in a static random potential corresponding to a height 
profile {hi}. The Hamiltonian is 

n = J2 h ^- (8) 

i 

Height differences are related to the hopping rates of the 
model through the condition of detailed balance 

exp[-(h i+1 - hi)/T] = (9) 

u i+l,i 



In this subsection we discuss both the steady state dy- 
namics and the coarsening dynamics of systems with uni- 
directional and bidirectional disorder. 

With unidirectional disorder, the behaviour of den- 
sity fluctuations in steady state resembles that in the 
zero-range process discussed in the previous section. For 
\p — p c \ > A, there is a kinematic wave which circulates 
through the system, and induces long-period oscillations 
in the density |Tfjj . The dissipation was monitored by 
performing a Galilean shift to keep up with the wave. As 
expected on the basis of the kinematic wave criterion dis- 
cussed earlier, disorder was found to be irrelevant, and 
z = I still holds ^(J- For p — |, on the other hand the 
system is phase separated, and the system supports two 
kinematic waves which move in oppsite directions in the 
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two coexisting phases. Simulations suggest that the kine- 
matic waves emanate from one interface, and terminate 
at another, leading to a different universality class (lfjj . 

In the phase-separated regime, the kinetics of approach 
to the steady state involves a coarsening process 0. 
Large-density stretches accumulate behind bottleneck re- 
gions which are stretches of consecutive weak bonds. If 
A J is the difference of currents exiting from two such 
bottlenecks separated by a distance £, the region in be- 
tween would fill up in time t ~ £/AJ. Estimating the 
probability of occurrence of stretches of length I and re- 
calling that they carry a current Ji — UijA^Y + a/l) 
leads to the estimate 

Turning to bidirectional disorder with tilt (/ < J), 
recall that in steady state, particles are held back be- 
hind the longest backbend, while the region ahead of it 
is empty. The coarsening properties of this model were 
obtained by Krug who pointed out that bottleneck 
stretch lengths are much smaller than their separations, 
implying that we may mentally replace stretches by sin- 
gle bonds that allow a maximum current Ji ~ A^' to flow 
through. Following through, we find that the coarsening 
length is given by [f| 

£(i)~£TW3. (12) 

In the untilted case, the steady state is the equilbrium 
state of hard core particles in a random (Sinai) poten- 
tial. Not much is known about the dynamics in this 
case, though the single particle problem is well studied 
|4lj . The existence of large (~ O(vT) potential barriers 
would make the approach to equilibrium logarithmically 
slow, and it would be interesting to see how much slower 
the coarsening process is than the familiar diffusion of a 
single particle in a Sinai landscape. 

IV. DYNAMICAL DISORDER 

In this section, we study the properties of an ASEP 
subjected to bidirectional bond disorder, where the dis- 
order is itself time-dependent. Specifically, we consider 
an untilted potential (as in Fig. 2(d)) which is evol ving 
stochastically in time through the single-step model |42| 
which preserves the overall no-tilt condition under time 
evolution. The basic dynamical move is the interchange 
of nearest neighbour links of the height field at a certain 
rate if they represent opposite slopes. If L and R stand 
for left-leading and right-leading links respectively, then 
a link configuration ....LR.... can evolve into ....RL.... at 
a certain rate, and the reverse move is allowed at an 
equal rate. The large space and time properties of the 
height field of this model are described by the continuum 
Edwards- Wilkinson model The ASEP particles are 
passive sliders on the surface, in that they are influenced 



by, but do not back-influence, the dynamics of the fluc- 
tuating surface. 

This problem has been studied in [3q|. Particles tend 
to form clusters when they fall into valleys, and to 
decluster when valleys overturn and form hills. How- 
ever, declustering involves activation and is a slow pro- 
cess, while clustering is relatively rapid. This asymmetry 
leads to overall clustering. Numerical simulations show 
that this tendency to cluster persists at all length scales, 
and leads finally to a steady state that exhibits phase sep- 
aration. However, this phase separated state is unusual 
in that is characterized by large fluctuations of the order 
parameter, which do not damp down in the thermody- 
namic limit. Despite the presence of strong fluctuations, 
the system never loses its ordered character |44|. Fluc- 
tuations carry the system relatively rapidly from one or- 
dered configuration to another macroscopically distinct 
one. But the probability for the system to leave this 
attractor of ordered states vanishes exponentially in the 
system size. 

A quantitative measure of the ordering is provided by 
the two point correlation function. Results of simula- 
tions show that Eq. 1101 derived for the disordered prob- 
lem with an untilted, bidirectional potential describes the 
form of Cl (x) in this case as well. When the update rates 
for the landscape and particles are equal, the value of 
m 2 is close to 0.8, indicating long-range order, consistent 
with the existence of a phase ordered state. Moreover, 
the value of the cusp exponent a is close to \ , its value 
for the problem discussed in the previous section. This 
quantitative similarity between the scaling properties of 
phase ordering induced by quenched disorder, and by dy- 
namical disorder corresponding to an evolving Edwards- 
Wilkinson landscape, is not completely understood. 



A. Time-dependent properties 

The time dependence of the particles is dictated by 
fluctuations of the driving height field. Except at 
very small times, the autocorrelation function A^it) = 
(rii(0)ni(t)} is found numerically 0] to be a scaling func- 
tion 

A L {t) = X(t/L z ) « m 2 [l - d{t/L z )% t/L z « 1 

(13) 

with z = 2 being the dynamical exponent governing the 
height-field dynamics and d a constant of order unity. 
Simulations show that m 2 ~ 0.8, /3 ~ 0.2. Like its coun- 
terpart for static correlations (Eq. I10[l . the scaling func- 
tion X shows an approach with a cusp singularity to a 
finite value at small argument, once again signalling long 
range order of an unusual sort. Finally, to understand the 
coarsening properties of this model, recall that cluster- 
ing occurs by collecting particles in valleys of the driving 
height field. Since the largest valleys that form in time t 
are typically of size t 1 ^, we expect the coarsening length 
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in this case to follow 

£(t)~t*. (14) 

The time-evolving two-point correlation function 
(ni{t)rii +x {t)) — p 2 is expected to be similar to the form 
of Eq. except that the system size L is replaced 

by the coarsening length £(t). The results of numerical 
simulations agree well with this [3^ • 

V. CONCLUSION 

We have reviewed the properties of the one- 
dimensional ASEP with various types of disorder - 
particle-wise, space-wise (unidirectional and bidirec- 
tional), and dynamical. All forms of disorder bring about 
phase separation, though under different conditions, and 
of different types. 

With particle-wise disorder, phase separation occurs 
below a critical density p c = (1 + /5 C ) _1 . The macroscopic 
stretch ahead of the slowest particle is totally devoid of 
particles (p = 0), whereas the stretch behind is mixed 
and has p — p c . 

With unidirectional space-wise disorder, the system 
supports phase separation in the region p c < p < 1 — p c , 
with coexisting phases that have densities p c and (1 — p c ). 

With bidirectional disorder with tilt, the system phase 
separates at all values of p into phases with densities 
close to and 1. The deviations from these values are 
proportional to the value of the current (Eq. |HJ), and 
approach zero as L — > oo. 

With bidirectional disorder and no tilt, we have the 
problem of hard core particles in a Sinai potential. The 
phases are again characterized by p = or 1 , as evidenced 
by m = 1 in Eq. I1UI but have some unusual properties: 



the interfacial region is very broad, and there are strong 
sample to sample variations in the length of the longest 
ordered stretch. 

Finally, the state reached through dynamical disorder 
is surprisingly similar to the previous case, with the dif- 
ference that m a =/= 1 in this case. 

Turning to dynamical properties in the steady state, 
we have seen that results are consistent with the kine- 
matic wave criterion |l6j . namely that a finite velocity 
of the kinematic wave renders disorder irrelevant; the 
universality class characterizing the decay of fluctuations 
remains unchanged. The determination of the disorder- 
determined universality class when the wave speed van- 
ishes remains an open problem, for instance, for particle- 
wise disorder with n < 1. 

The coarsening behaviour characterizing the approach 
to the steady state is well understood in several cases 
(Eqs. 131 1111 IT21 114|) . However, the coarsening properties 
of a system of hard core particles in the Sinai potential, 
which corresponds to the untilted bidirectional case, re- 
mains an interesting open problem. 
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